Analytical progress on symmetric geometric discord: Measurement-based upper 

bounds 
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Quantum correlations may be measured by means of the distance of the state to the subclass of 
states n having well defined classical properties. In particular, a geometric measure of asymmetric 
discord [Dakic et al., Phys. Rev. Lett. 105, 190502 (2010)] was recently defined as the Hilbert- 
Schmidt distance of a given two-qubit state to the closest classical-quantum (CQ) correlated state. 
We analyze a geometric measure of symmetric discord defined as the Hilbert-Schmidt distance of 
a given state to the closest classical-classical (CC) correlated state. The optimal member of is 
just specially measured original state both for the CQ and CC discords. This implies that this 
measure is equal to quantum deficit of post-measurement purity. We discuss some general relations 
between the CC discords and explain why an analytical formula for the CC discord, contrary to the 
CQ discord, can hardly be found even for a general two-qubit state. Instead of such exact formula, 
we find simple analytical measurement-based upper bounds for the CC discord which, as we show, 
are very efficient in the case of two qubits and may serve as independent indicators of two-party 
quantum correlations. In particular, we propose an adaptive upper bound, which corresponds to the 
optimal states induced by single-party measurements: optimal measurement on one of the parties 
determines an optimal measurement on the other party. We discuss how to refine the adaptive 
upper bound by nonoptimal single-party measurements and by an iterative procedure which usually 
rapidly converges to the CC discord. We also raise the question of optimality of the symmetric 
measurements realising the CC discord on symmetric states, and give partial answer for the qubit 
case. 



I. INTRODUCTION 



Entanglement is a fundamental type of quantum corre- 
lation that has come to be seen as an important resource 
in Quantum Information (see, e.g., Ref. [ij). However, 
quantum mechanics supports other, distinct from entan- 
glement, types of quantum correlations in composite sys- 
tems, such as the so-called quantum discord [JIH, whose 
characterization is the topic of much current research (see 
the review and references therein). Quantum discord 
is an information-theoretic measure of correlations where 
quantum correlations are identified in terms of the dif- 
ference of two classically equivalent definitions of mutual 
information [1, Q in a composite system. A different 
possible perspective on quantumness of correlations is 
captured in terms of quantum deficit functions Q, i.e., 
differences between certain properties of a state, before 
and after classical type measurements are performed on 
it. One such important property is the optimal ther- 
modynamic work that can be extracted from a state in 
scenarios of classical (local) measurement complemented 
by zero-, one- and two-way classical communication be- 
tween measuring parties Q (a state is classical if the 
deficit is zero). While the two-way scenario is rather in- 
volved, the zero- and one-way quantum work deficits are 
simply equal to the so-called relative entropy of quan- 
tumness jB) [3 - the minimal entropic "distance" measure 
to specific classes of classical-type states. 

Distance measures to sets of states with only classi- 



cal correlations are promising, and conceptually, simple 
ways of identifying quantum correlations. Recently, e.g., 
Dakic et al. [8] introduced a geometric measure of dis- 
cord of a state as its minimal Hilbert-Schmidt distance 
metric to the set of states with null quantum discord 
(these states are one-side classical, or so-called classical- 
quantum (CQ) states of the form p — '^.^PiPi® Pi, where 
Pi's are orthogonal projections with rank one and piS are 
quantum states). 

A natural, symmetric measure of quantum correlations 
can be obtained by constraining to a set of fully classical 
states, i.e., classical-classical (CC) states which are di- 
agonal in some product basis [9|. The optimization pro- 
cess required in the evaluation of (general) quantum cor- 
relation measures renders their calculation challenging. 
Here, we shall build on an equivalence between geomet- 
ric measures of quantum discord and quantum deficits of 
purity to provide tight and faithful upper bounds on the 
symmetric geometric discord. 

The paper is organized as follows. In Sec. II, we pro- 
vide some basic definitions and theorems for the discords 
in relation to quantum deficit. In Sec. HI, we present our 
main result - the measurement-based upper bounds on 
the CC discord. In Sec. IV, we give explicitly formulas 
for the upper bounds in the case of two qubits. In Sec. 
V, we present an analytical comparison of the discords 
and upper bounds for some classes of states. In Sec. VI, 
we present a few methods with examples for optimization 
of the upper bounds. We conclude in Sec. VII. 
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FIG. 1: (Color online) Venn-type diagram showing the sets of 
the CC {Q.S = ^ab), CQ (Qa), and QC {Q.b) states together 
with the closest states a* (i = S, A, B) according to the CC 
[Ds), CQ {Da), and QC {Db) geometric discords, respec- 
tively. States as' and (Tg//, which correspond to the adaptive 
upper bound Dg'"''' , are the closest CC states for cj\ and o*g , 
respectively. This is an intuitive graph but, more precisely, 
the point a*^ (cts) should be on the line between p and agi 
(c^S"). 



II. BACKGROUND 

We start by recalling the quantum zero-way and one- 
way work deficit 6]. Let the sets of states ft a, ^b, 
and fls correspond to classical-quantum (CQ), quantum- 
classical (QC) and classical-classical (CC) states, respec- 
tively (see Fig. 1). Note that the set fis = ^ab — 
n J7b is obviously in the intersection of the other two 
sets, and any element of the intersection is in the set. 
Let Mx correspond to all von Neumann's measurements 
that are associated with the set fix {X = A, B, S) in the 
following natural way, Ma —Ma®Ib,Mb ^ Ia®M.b 
and Ms = Ma ® Mb, where Ma, Mb are just local 
von Neumann's measurements performed by Alice and 
Bob in some orthonormal basis. We define the corre- 
sponding one-way {X = A oi B) and zero- way {X — S) 
quantum work deficits as: 

^x{p)^rn\nS[Mx{p)]-S{p), (1) 

Mx 

where S{-) is the von Neumann entropy. The relative 
entropy of quantumness 0| is 

D^{p) ^ uAn S{p\\a). (2) 

There is an observation (see Sec. VLD in Ref. @) that 
links the above quantities. 

Observation 1 . — For any quantum state it holds that 
[p) = [p) . Furthermore 

D'^{p)=uAuS{p\\Mx{p)), (3) 

Mx 



which means that the optimal state a*x, saturating the 
minimum in Eq. (0), comes from the optimal measure- 
ment in Eq. ^j) of the examined state M*x{p) — cr*^- A 
proof of Observation 1 is given in Appendix [X] 

Lemma 1. — For any function f , any Hermitian opera- 
tors F and G, and any von Neumann's measurement op- 
eration M we have tv[Ff{M{G))] = tr[X(F)/(X(G))]. 
See Appendix |X] for a proof of this Lemma. 

Geometric discord as a purity deficit. — Geometric 
measures of discord of a state are also similarly com- 
pletely determined by optimal measurements on it, as 
already noted by Luo and Fu [l^l and elucidated in the 
review by Modi et al. . We formulate this property as 
follows: 

Observation 2.- Let a*^ G VLx he an optimal state sat- 
urating the minimum for quantum geometric discord [^]: 

Dx{p)^ mm\\p-a\\' (4) 

defined by the norm \\A\\ — ^/tT{A^A). Then it is re- 
alised by some optimal measurement M*x on p, i.e., 
a*x — M*x{p) and satisfies the Pythagorean formula: 

Dx{p) = \\p-o*x\? = \\pf-\\a*x\?. (5) 

Thus, M*x maximizes the post-measurement purity 
max^^ iT:[{Mx{p)Y] leading to the alternative formula: 

Dx{p) = vainWp- Mx{p)\\^ 

Mx 

- tr(p2)_maxtr[(^x(p))2]. (6) 

Mx 

See Appendix El for a proof of this observation. Note 
that choosing X ^ A corresponds to one-side or asym- 
metric (CQ) geometric discord Q, while X = S corre- 
sponds to the symmetric (CC) version Q. The last form 
Eq. ^ for geometric discord highlights an immediate 
analogy to the original deficit of Eq. ([IJ on replacing the 
original von Neumann entropy S{p) = Sa=i{p) by the 
Tsallis entropy Sa{p) = --^tT:{p°') (for a = 2). We 
shall refer to the left-hand-side of Eq. ([6]) as a purity 
deficit, which is a special case of entropy based deficits 
AJ'^ :— min_^^ Sa[Mx{p)] — Sa{p)- For completeness, 
we provide a proof of Observation 2 in Appendix |^ (see 
Ref. [131 and Sees. ILG and IILB.2 in Ref. 3 for aher- 
nate proofs). 

Simple consequences. — Observation 2 leads us to: 
Lemma 2.- For both the geometric discords Dx and 
quantum discords Dx based on relative entropy, it holds 
that: (i) The CQ and QC discords bound the CC dis- 
cord from below, Ds{p) > raax.[DA{p), DbIp)]- (H) We 
have simple implications Da{p) = Ds{p) = Db{p), 
and Db{p) = =^ Ds{p) = Da{p) with optimal mea- 
surement M*g being a product of the measurement re- 
alising the respective CQ or QC discord with the one 
that commutes with the initially classical subsystem. (Hi) 
Ds{p) = ^ Da{p) = Db{p) = 0. 

See Appendix [X] for a proof of this Lemma. 
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III. MEASUREMENT-BASED UPPER BOUNDS 
ON THE CC DISCORD 

We now turn to the main result of this paper. The ex- 
pUcit calculation of the geometric CC discord is in general 
difficult as it involves optimization over all measurements 
of the required form given by Eq. ([6]). In particular, the 
CC discord involves optimization over two sides of the 
states and so involves twice as many parameters as the 
CQ case. For the case of CQ-type discord, certain lower 
bounds have been found ^Oj, JJj . On the other hand, we 
show here that the measurement based formula Eq. ([6]) 
can be fruitfully used to construct useful upper bounds 
on the CC discord. 

The conclusions of the present paragraph are valid for 
the geometric and relative entropy discords and for quan- 
tum deficit based on any quantum entropy Sa- 

Let us recall that one refers to a bound as: (i) tight if it 
coincides with the bounded quantity on some non-trivial 
subclass of states and (ii) faithful iff it vanishes on any 
state for which the bounded quantity vanishes. 



A. Nonadaptive upper bound 

An arbitrary measurement over two sides of the state 
is, by definition, an upper bound on discords: 

jj^=Ms(P)^p^ = S^ia)-S^{p) (7) 
= S^[Ms{p)]-So.{p)>Ds{p), 

where a G [0, oo]. For ease of notation, let Ai*x p denote 
the optimal measurement leading to the discord Dx of 
state p. Product of the two (CQ and QC) optimal mea- 
surements on the state p leads to the first interesting 
bound, that we shall call the simple product ( or nonadap- 
tive) bound, for which the measurement-induced state is: 

^=[MXp®M*bJ{p) (8) 

in Eq. ([7]). This is one of the simplest kinds of bounds 
motivated by asking how the CC and CQ discords (or 
optimal measurements) are related. Indeed, we have al- 
ready noted in Lemma 2 that this type of bound trivially 
coincides with the CC discord in the special case when 
one of the CQ discords is null. 



Note that part (ii) of Lemma 2 immediately leads to the 
following 

Fact 1. — The hounds ^ based on measurements, 
given by Eqs. 0), and HU\) . are faithful, so they 

may serve as independent indicators of two-side quantum 
correlations. 



C. Iterative procedure for the adaptive upper 
bound 

The adaptive form of Eqs. ^ and (|10p allows for an 
iterative procedure that may be helpful in refining up- 
per bound on the CC discord. Indeed, let X and X' 
be two opposite subsystems [i.e., {X,X') = {A,B) or 
{X,X') = {B,A)]. Consider the following procedure: 
Step 1. Choose the initial subsystem X = Xq (cither A 
or B), and initial measurement Aix — ■^*Xo- ^^^P ^• 
Iterate the following steps: Step 2.1. Given input mea- 
surement Mx on X calculate the output, i.e., optimal 
measurement Ai*xi Mx(p) second system X' . Step 

2.2. Put X' in place of X and the output M*x, Mx(p) 
the input for Step 2.1, calculate its output again. Step 

2.3. Calculate the bound on the discord, given by Eq. ([7|), 
with the help of the measurement AAs being the tensor 
product of the input-output pairs of measurements on X 
and X' presented in Steps 2.1 and 2.2, take the minimum 
of the two. Step 2.4- Take the minimum of the output 
of Step 2.3 of two subsequent rounds. Step 3. Stop the 
procedure if the outcome of Step 2.4 does not change. 

IV. TWO-QUBIT CASE REVISITED 

We now consider the case of the CC discord of two 
qubit states. The standard Bloch representation of any 
two-qubit state is 

1 ^ 
p ^ -{I(g)I-\-x-a(g>I-\-I(g>y-a+^ T^ja^ (g> aj) 

= f{x,y,T) = f{\x),\y),T), (11) 

where a — [01,02,0^] is a vector of three Pauli matrices, 
T is the correlation matrix with elements Tij = tv[p{ai ® 
CTj)]; f = [xi,x2,X'i]^ = \x) and y = [yi, y2, 2/3]^ = \y) 
are the (column) local Bloch vectors with components 
Xi = iv[p{ai (g) 1)] and yi — tr[p(Il ® Ui)]. 



B. Adaptive upper bound 

One can further introduce refined bounds that are 
adaptive, i.e., measurement on one of the parties is per- 
formed on the optimal state corresponding to the other 
party, as below: 

^" -[MXm',Jp)®M*sJ{p). (10) 



A. CC vs CQ discords 

We state the following simple 

Fact 2. — Any two-qubit state p — f{\x),\y),T) is 
mapped into 

^in)Ap) = f{\h){n\x),\y),\h){h\T), (12) 
^(a)b(/') = f{\x),\i^{My)-,T\m){m\), (13) 
(^(n,rh)s{p) = /(|n)(n|a^),|m)(m|y),|n)(n|r|TO)(TO|)(14) 
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by the measurement of (i) ha on the left quhits, (ii) rha 
on the right qubits, and (in) ha and rha on the left and 
right qubits, respectively. 

This follows from Lemma 1 and the fact that diagonal 
of ha vanishes in the eigenbasis of any h'a with n' _L n. 

Observation 2 and Fact 2 directly lead to the analytical 
formula (see Ref. [8|) for the CQ discord Da as follows: 



DaIp) 



l|2 



tr(p2) _ niax{tr[(a(^) Jp))^]} 

k 



\x\\^ + \\T\\^ - ma^[{k{\x){x\+TT^)\k)] 



\T\ 



(15) 



where kx is the largest eigenvalue of matrix — |a;) 
TT^. For clarity, we also write 



4||p||2= l + P||2 + ||y||2 + ||T||2 

= l + {x\x) + {y\y)+tT{TTT). 



(16) 



However Observation 2 yields more, viz. the eigenvector 
\kx) corresponding to the eigenvalue k^ defines the opti- 
mal measurement of party A on p producing the closest 
CQ state a, which {via the Fact 2) is 



c^l = f{{kx\x)\K), \y), \kx){kx\T). 



(17) 



Analogously, one obtains Db{p) = \{tT:Ky — ky) — 
idlyiP + II^IP ~ ^y)j where ky is the largest eigenvalue 
of matrix Ky = \y){y\ + T'^T with the eigenvector \ky). 
The closest QC state is 

a*s= f{\x),{ky\y)\ky),T\ky){ky\). (18) 

Observation 2 also delivers the two-qubit CC discord 

^5(p) = ||p|p-|kSlP (19) 

with the norm of a*g = /(|x^), \y*s),Tg) can be given in 
terms of some functions minimized solely over unit vec- 
tors \xs) (or, equivalently, \ys)) b& given in Appendix IB] 
This is identical to the single Bloch-sphere optima ob- 
tained in Ref. [1]. 



B. Quest for symmetry of the optimal 
measurement for symmetric states 

There is a general question whether the states symmet- 
ric under swapping subsystems always allow for a sym- 
metric optimal measurement in the formula for the CC 
discord Ds- Here, we provide some partial results on this 
problem. Namely, there is a practical observation: 

Theorem 1 Consider the two-qubit symmetric states p, 
i.e., the ones satisfying pAB = PBA or, equivalently, 



T = 
\x) 



\y)- 



(20) 
(21) 



// the matrix T satisfying either T > or (— T) > then 
the optimal CC state crj and the corresponding measure- 
ment are .symmetric, i.e., the optimal measurement basis 
is defined by some \x*g) — \yg). 

Proof. — Clearly, since Ds{p) = tr(p^) — 
maxig^yg Tr[cr(£g.j;g)^j5 (p)^], we may write it in the 
form 



7^5(p)=tr(p2)_i 



1 + maxM(x5,ys) 

XS,VS 



where the function u is defined as 



u{xs, ys) = {xs\T\ys)'^ + {xs\x)'^ + {vs\vf 



(22) 



(23) 



Following Theorem 1, it is not difficult to see that, by the 
symmetry of the initial state p, one has |a;) = \y). Now 
for T > (all eigenvalues strictly positive) one defines 
the new scalar product {xs,ys)T — {VTx\\/Ty) , which 
defines also the norm ||ss||t — \/ [xs, xs)j.. Now, since 
\\xs — VsW't — 0, for any pair of unit vectors \xs) and 
\ys) one has ^[u{xs,xs) -{-u{ys,ys)] > u{S:s,ys), which 



means that the maximum in Eq. (|22p is achieved by a 
symmetric pair {\x*g) = lyj)). The proof for (— T) > 
goes along the same lines. For the cases when T > 
or (— T) > 0, i.e., where zero eigenvalues are allowed, 
the statement follows from the continuity argument since 
here the argument realizing the maximum is continuous 
in parameters of the state. 

We conjecture that in the case of the symmetric two- 
qubit states any minimum can be reached by symmetric 
measurement. We have performed both analytical and 
numerical search and found no counterexample to this 
hypothesis so far. However for higher dimensions it may 
not be true since as we know there are numerous prop- 
erties that break there. 



C. Adaptive and nonadaptive upper bounds 

We now turn to the upper bound for the CC dis- 
cord. Using the adaptively measured states, given by 
Eqs. (IH) and (fTO|) . we obtain the following upper bound 
from Eq. Q for a = 2. 

Theorem. — For an arbitrary two-qubit state the adap- 
tive upper bound D^g^^\p) for the CC discord can be 



given by: 



D 



(aub) 



ip) 



mm p-CTi = 



- max||cr,;|p,(24) 



where the CC states agi and as" are 



as' = f{\xs'),\ys'),Ts')^a^^,^i^^^{p) (25) 

- f{{kx\x)\k)Aly\y)\k)Akx){k\T\ly){ly\), 

^S" = f{\xs")Ays").Ts'') = a^^i^f^^^^{p) (26) 

= f{{ix\x)\ix)AK\y)\k)AQ{l\T\ky){ky\)- 
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where \k.j;) , \ky), \lx), and \ly) are the eigenvectors corre- 
sponding to the maximum eigenvalue of 



Kx 



xx^ +TT^ EE |a;)(a;| 
\y){y\+T^T, 



= \x){x\+T\ky){ky\T'^, 

Ly = \v){v\+T'^\k){K\T. 



(27) 
(28) 
(29) 
(30) 



respectively. Note that as' in general differs from a*g 
used in Eq. p9p . Explicitly, the norms are given by 

11^5'!^= \{i + \\xs'\? + \\ys'\? + \\Ts'\?) 

^ \{l + {kW + {iy\yf + {k\T\ly)^), (31) 
1 1^5" IP = + + \\ys"\? + \\Ts"\?) 

= \{l + {l\x)^ + {ky\y)^ + {l\T\ky)^), (32) 

Note that the measurement on direction \lx{y)) corre- 
sponds to the adaptive measurement Af^j-^-j ^, since the 
correlation matrix of the optimally measured state on 
A {B) is according to Eq. dTT]) [Eq. ([18])] given by 
\kx)(kx\T (r|fcj,)(fcj,|). Intuitively, the adaptive upper 

bound D^s"''\p) can also be found by applying the fol- 
lowing relation between the three discords valid for an 
arbitrary two-qubit state: If Da{p) — then Ds{p) = 
Db{p) and, analogously, if D_b(p) — then Ds{p) = 
Da{p) as given by Lemma 2. The bound can be con- 
structed as follows (see Fig. 1): 



4^"^'(p) = min(i?s',i?S"), 



(33) 



where 

Ds' ^\\p-<J*a\ 
Ds" = \\p-cr*B\ 



Iks" 



(34) 



and as' and as" , given by Eqs. (1^ and (1^ . were calcu- 
lated from the repeated application of Eqs. ([T7| and ([T8| . 
It is also worth noting that 

^s"''^(p) = ^ Dsip) = ^ Da{p) = Db{p) = 0. 

(35) 

So, in particular, it means thatZ?^""'^)(p) is nonzero iff p 
is not a CC state, and thus it may serve as an indicator 
of quantum correlations itself. 

The nonadaptive upper bound (i.e., product bound) 
for a two-qubit state p can be given by 

D^r''\p)^\\p\?~\\asA\\ (36) 

where 

'TSo =/((fcxk)|fc.),(fcyl2/)|fcy),|fcx)(fcx|T|fc,)(fc,|), (37) 
for which 

IksolP = ^(1 + {k\xf + {ky\yf + {kx\T\kyf). (38) 



We have the following inequalities 



max(i^^,I?B) <Ds< D^;^""^ < 



(39) 



where the last inequality can be immediately concluded 
by comparing Eqs. ^ and with Eq. ((M)) . 

We note here that the adaptive bound, given by 
Eq. ([M)) . is very effective. Indeed, the largest gap to 

the exact value S — D'g^^\p) — Ds{p), observed by us 
numerically, is just a few percent, and it is usually of the 
order 10^^ or 10~^ for randomly generated rank-4 states. 
Interestingly, we have also observed that it is exactly zero 
for almost all classes of states for which there are known 
analytical expressions for Ds- 



V. DISCORDS AND UPPER BOUNDS FOR 
SOME CLASSES OF STATES 

A. Examples of simple relation between discords 
and their upper bounds 

Here, we present some examples of analytical calcula- 
tion of the CQ and CC discords and the adaptive upper 
bound to show their relations. 

Example 1. — For (a) pure states, (b) Bell diagonal 
states, and also for (c) states with both marginals van- 
ishing, i.e., \x) = \y) = 0, it holds 



D 



(aub) 



(40) 



For these states, Ds can be easily found by showing ex- 
plicitly that the lower bound Da = -D_b is equal to the 
upper bound D^g^^K 

Example 2. — For states with maximally mixed single 
marginal, e.g., \x) = (and analogously for \y) = 0), we 
have 



4||cr^ 



= 1 



{x\x*sy +(y\y*s) 

=0 



{is\T\ys) 



Since \x*g) maximizes ||tTg 



l + {y\ysr + {i*s\[T\ys){ys\T^] Ws)m 
P then it holds 

T\rs) 



ViyslT^m) 



(42) 



Thus, we obtain 



{ys\T' 



T\y*s) 



1 + {y\y*sf + {ys\T^T\y*s) 
i + {y*s\{\y){y\+T^T)m 

l + max[eig(|y)(yH-T^T)], 



(43) 
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FIG. 2: (Color online) Geometric discords and their tight 
upper bounds for the state p{p, (j> ~ 7r/2), given by Eq. (|46|) . as 
a function of parameter p: The CC discord, Ds (blue solid), 
and CQ/QC discords. Da = Ds (magenta dotted curve). 



together with the adaptive upper bound, D^^"''' (red dashed), 



(nub) 



and unoptimized nonadaptive upper bounds, i'gn 
(green dot-dashed). We note that D^g^^^ 



D 



(nub) 
S22 



D 



(nub) 



p,{nub) 



Ds and L»lf"''^ 



D 



321 



for any p. All the upper bounds, except 



S12 

P){aub) 
S 

vertical lines are added for clarity only. 



D''a'^"' , are discontinuous at p = 1/2, while the corresponding 



SO finally 



Ds 



^{{y\y) 



\T\ 



max 



[eigi\y){y\+T^T)]} 



(44) 

which is equal to the QC discord Db and the adaptive up- 
per bound D^g^'"^ , which follows from a simple direct cal- 
culation. By performing analogous derivation for \y) = 0, 
we conclude that 



la;) =O^Da<Db = Ds = 

\y) =O^Db<Da^Ds^ D^;''''\ 



(45) 



B. 



Example of nontrivial relation between discords 
and their upper bounds 



Here, we give an example of nontrivial relation between 
CC, CQ, and QC discords and their tight upper bounds 
as shown in Figs. 2 and 3. Specifically, we will study 
mixtures of BeU's state l^-^) = [|01) -I- exp(i(/))|10)]/v^ 
and |00) (i.e^ state separable and orthogonal to \'^^)) as 
defined by [l|,[l3: 



p(p,0)=p|*^)(*^| + (l-p)|OO)(OO| 



(46) 



FIG. 3: (Color online) Same as in Fig. 2, but for the state 
p{p, (j)) as a function of phase (j) for fixed p — 2/3. 



for < p < 1. Wc find that the CC discord for these 
states is given by 



Ds = - min [2p^ , 7p^ 



P + 3] 

if P<f, 

(7p^ — 8p + 3) otherwise. 



(47) 



By contrast, the CQ and QC discords are given by 
1 



Da^Db = - min(p^ , - 3p + 1) 



if 



3p+l) 



P 



< 



otherwise. ^^^-^ 



Some details of the calculation of the discords are given 
in Appendix [C] Moreover, we find the adaptive upper 
bound for the state p{p, 0) to be 



D 



(aub) 



Ds'=Ds" = 



(aub) 



3) otherwise.^ ^ 



Ds for < p < 1/2 and 3/5 < 



It is seen that D 
p<l. 

The nonadaptive and adaptive upper bounds for this 
state can be optimized as will be described in the next 
section. All these discords and upper bounds are shown 
in Fig. 2. In particular, we observe discontinuity of the 
upper bounds at p = 1/2. We find that the upper bound 



D^g^^^ (and I?^"""^) has two different limits: 



,(nub) 



hm D 



(aub) 



1 

8' 



Um D 



(aub) 
S 



3 

16' 



(50) 



By contrast, the asymmetric discords Da — Db, sym- 
metric discord Ds and optimized upper bounds D^^^^^ 
(as discussed in the next subsection) are continuous func- 
tions of any p. Anyway, none of the discords has contin- 
uous first derivative in p. 
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VI. IMPROVED UPPER BOUNDS 

A. Optimization over degenerate measurement 
outcomes 

If the maximal eigenvalues of operators Kx,y and/or 
Lx^y are degenerate then the adaptive and nonadaptive 
upper bounds can be optimized by taking the minimum 
for the eigenvectors corresponding these maximum eigen- 
values. Here, we will describe this method in brief and 
give an example explaining Figs. 2 and 3. 

First, it is worth recalling now a classic linear- algebraic 
theorem stating that eigenvectors of degenerate matrices 
are not necessarily orthogonal, but they can be made 
orthogonal and complete, as in nondegenerate case, by 
applying Gram-Schmidt's orthogonalization procedure. 
This is possible by having additional freedom of replacing 
the eigenvectors corresponding to a degenerate eigenvalue 
by their linear combinations. 

Let us denote eigenvectors \kx^) {\ky^)) corresponding 
to the same maximum degenerate eigenvalue of operator 
Kx {Ky), given by Eq. ^ [Eq. (gS])]. Analogously, we 

denote eigenvectors \ix''^) and \iy''^) corresponding to the 
maximum degenerate eigenvalues of operators: 

L« = |x)(x|+T|^W)(^W|rT, (51) 
L« = |z/)(z/|+T^|fcW)(fc«|T, (52) 

respectively. Thus, by applying these eigenvectors to 
Eqs. (ISIl), and ^E^, one can obtain norms ||o-^'|;;'''|P, 
||(Tg/'|p and ||o-5//''|P resulting in: 

^st^ - „min„(llp|P-|k^'^'^in, (53) 
= IHP-lkg^ir (54) 

Then, the optimized adaptive and nonadaptive upper 
bounds are simply given by 

Z?'""*^) = min , Z^f""*^) = min Z^^'jf ^ , (55) 

respectively. 

Example. — Let us analyze again the state p{p,4>), 
given by Eq. PS|) . Operator Kx is degenerate, as given 
by Eq. (|C2[) . so we can choose \kx^) = \i). Simple cal- 
culation shows that one can also choose ifi*"''') = |j) for 
i,j = 1,2. We find that the nonadaptive upper bound 
for z = 1 , 2 is equal to 

^(nub) _ f f if P < 5, (.r.. 

\ |[p2(cos2,/.)-8p + 3] otherwise. ^ ^ 

as shown by the green curves in Figs. 2 and 3. By con- 
trast, i?^^^''^ = Df^^^ = Z?^""^) as given by Eq. (03. 
So, finally, 

D^""^) = min(i?^™'^\i?^'i"2^') = d'^^^^ = i?^""^'(57) 



as shown by the red curve in Fig. 2. Note that such 
degenerate-value optimization for D^g^^^ is unnecessary 
for this state. 

By analyzing our formulas and Fig. 2, we can observe 
that (1) Da = Db ^ Ds for p € (i,l), (2) dIj") ^ 
for p e (i, 1) if ^ 0, and (3) i?^""^^ = D^p""^ = 

O'g^^^cj) = 0) 7^ L»s for p e (i,f). We observe that 
the unoptimized nonadaptive bound can be much greater 
than the adaptive bound if ^ ^ and i < p < 1, thus 
including the case for Bell's states {p = 1). In Fig. 3, we 
analyze the state p{p,(f>) for p = 2/3. We observe here 
that (i) the symmetric discord (blue solid line) is equal to 
the adaptive upper bound, Ds = D^g^^^ — 7/36, (ii) the 
asymmetric discords (black dot-dashed line) are Da = 
Db = 1/6, (iii) the nonadaptive upper bound (red dotted 

curve) depends on (j) as follows D^s^^^ = Ds + sm^{(t>)/^. 

Finally, we conclude Ds = i^^"""^^ = D^r''^ = D^^^'^ < 

D^s^^'' for < p < 1/2 and 3/5 < p < 1. We see that the 
nonadaptive bounds without optimization, on the other 
hand, can fare rather badly as an estimator of the CC 
discord. 



B. Optimization by locally nonoptimal 
measurements 

Here, we suggest to optimize the adaptive upper bound 
by locally nonoptimal measurements, i.e., to optimize 
over all measurement outcomes corresponding to all (for 
i,j ~ 1, 2, 3) measurements of party A (B) on p produc- 
ing (usually not the closest) state tr^^ (o'b'') and then 
measurements of party B (A) on this state producing the 

state agP (c^??). Thus, we describe the optimization 
of the adaptive upper bound over all eigenvectors of Km 
and L,„ (m = x, y) corresponding to all eigenvalues in- 
stead of taking only those eigenvectors corresponding to 
the maximum eigenvalues of these operators. This some- 
how counterintuitive method can be in fact efficient for 
the adaptive upper bound since L„i are constructed via 
eigenvectors of Km ■ Clearly, the method cannot improve 
the nonadaptive upper bound, as the operators Lm are 
not used there. 

The optimized adaptive upper bound Z)^™'^-' can be 
defined in analogy to Eq. ([55)1 as follows: 

min Ds,,^\\p\f~ max max \\ai'^^\f. 

i,i=l,2,3 •^'■^ " " r=S',S" i,j=l,2-3 

(58) 

By contrast to Eq. the optimalization is over 2x9 

parameters for any state independent of its degeneracy. 
It is convenient to form 3x3 matrices with elements 
will be done in the following. 
Example. — Again we analyze the state, given by Eq. 
For each of the three eigenvectors \km) of Km (for 
TO — x,y), given by Eq. \C2\ . we find three orthogo- 
nal eigenvectors \im ''), according to Eqs. ((?T|) and ([5^ . 
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Then, we can calculate ||(t 
e.g., the following matrices: 



112 
S' II 



" 1 

4 


A 


c 


1 

4 


A 


c 


c 


C 


B 


' c 


c 


B 



and create, 



(59) 



C A 
C A 



if 



(60) 



2x1 



where A = {1 + Tl^)/A = (1 + p2)/4, B = (1 
T|3)/4 = [1 + 2(1 - pf + (1 - 2p)2]/4, C = (1 + xl)/A = 
[1 + (1 — p)^]/4. Any order of the eigenvectors (and, thus, 
the order of the elements in the above matrices) can be 
applied. For convenience, we ordered them here by the 
value of the corresponding eigenvalues. Then, we obtain 



D 



(aub) 



-max(A, B, C) = 



B 
A 



< 



if P^5, 
otherwise, 



where \\p\\^ = '^p{p — 1) + 1 (see Appendix C). Thus, we 
conclude that the optimized upper bound is equal to the 
CC discord for any p £ [0, 1]: 



D 



(aub) 



= - min [2p^, 7p^ — 



in agreement with Eq. (|47| . Note that A 
and A = C for p = 1/2. So 

1/2 contrary to discontinuous D^""*^^ 



3] = Ds. 

B for p 



(61) 

= 3/5 



^^aub) continuous at 



P 



and D 



(aub) 



(compare broken and solid curves in Fig. 2). 

In conclusion, for p(p, ip) with l/2<p<3/5 and any 
(j), we have the following inequalities 



Da = Db <Ds = D 



(aub) 



< D 



(aub) 



= D 



(nub) 



(62) 



This example demonstrates usefulness of the optimiza- 
tion procedure by calculating the upper bounds for all 
possible measurements rather than only for those mea- 
surements corresponding to the maximum eigenvalues of 
Km and Lm {m ^ x,y). 



TABLE I: Examples of the application of the iteration pro- 
cedure for the adaptive upper bound D^g^^^ for the states 
given by Eqs. dSSI, ([TTJ, and as described in Sec. lVICI 

The accuracy of the procedure is shown by the difference A„ 
between the adaptive upper bound after the nth iteration. 



D 



(aub n) 



, and the exact value of the CC discord, Ds- 



iteration 


A„ 


T-^(aub n) 

— Ug - 


Ds 


n 


state (f68l) 


state (f7T|) 


state (1721) 





8.85x10-* 


4.28x10"^ 


1.71x10"" 


1 





4.77x10"^ 


8.00x10"'' 


2 




5.53x10"^ 


3.90x10"'^ 


3 




6.44x10"" 


1.92x10"* 


4 




10-" 


9.50x10"" 


5 




10-" 


4.69x10-" 



where \kx 

{n} 

y 



:{o} 



while 



and \ky) are the eigenvectors correspond- 



ing to the maximum eigenvalues of 

4"> = |x)(x|-f-r|fcM)(fcM|TT, 
= |y)(y|+r^|fci">)(fcM|T, 
ifW ^ K, = \x){x\+TT'^, 



y 



\y){y\+T^T, 



(64) 
(65) 
(66) 
(67) 



respectively. For randomly generated rank-4 states (thus, 

usually, with nondegenerate eigenvalues of K^^y and 

Li^y), the procedure is usually effective as can be shown 
be calculating the difference 

p.(aubn) 



between the adaptive upper bound after the nth iteration 
and the exact value of the CC discord. 

Let us discuss just a few examples: Example 1. — Let 
us analyze state p = f{\x), \y),T) described by: 



T-idiag([l, 



-1,01 



(68) 



Iterative procedure for the adaptive upper 
bound 



Here, we describe in detail the iterative procedure de- 
scribed in Sec. IILC for the adaptive upper bound Z3(^"^) 
and give some examples. The nth iteration of the adap- 
tive upper bound, D^^' 



.)(aubn) Qg^j^ l-,g (,j^J(,^lg^^g(J g^g 



D 



(aub n) 



where our old Dg^"'^'' is just 



^S' II ' W^S" II I 
(aub 0) 



(63) 



and 



First, we calculate the closest CQ state to be given in 
Bloch's representation as cr^ = /(|x^), \x),Ta), where 





' 1 + V3' 




xa) = t 


1 + V3 


, Ta = 




2 





t 
t 

1 

3\/3 



t 
t 




1 

8\/3 



(69) 



r{"}||2 _ 

^S' II — 

{"}||2 _ 



i(i-h(fci">ix)2 + ([riy)2-K(fcr^iT|[r)2) 

i(l -t- (f|">|x)2 + (fc^">|y)2 + (/l">|r|fc^">)2), 4""*^) = 5(3 - V3)/192 = 0.033- 



r,{"}|rTn|;{n}\2\ 



with t = {3 + \^)/A8. Analogously, the closest QC state 
is = f{\x), |xyi),rj). Thus, the CQ and QC discords 
are given by Da = Dg ^ {3 - V3)/64 = 0.0198 ••• . 
By contrast, the closest CC state is much simpler as 
given by a*g = f{\x),\x),Z), where Z is the zero ma- 
trix. Thus, the CC discord is simply equal to -D5 = 
1/32 — 0.031 ••• . The nonadaptive upper bound is 

which is obtained 
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as the Hilbert-Schmidt distance of p to the CC state 
t^So = f{\xA)AxA),Z), where \xa) is given in Eq. 
The CC states, defined by Eqs. (|25l) and (|26|) . are equal 
to as' = f{\xA)Ax),Z) and cts" = /(|a;), 2"), re- 
spectively. Thus, the (initial) adaptive upper bound is 
^(aub) _ ^o^ubo) ^ - 573)7384 = 0.032- ••. Our 
iteration procedure of the adaptive upper bounds con- 
verges to Ds already in the first iteration as Z?^^"*^ = 

1/32 (see Table I) since the CC state a^g} = a^^} = 
f{\x),\x),Z) = as- Thus, for the analyzed state, we 
have the following inequalities 

Da = Db<Ds^ D^s""" < < • (70) 

Example 2. — Another state is given by the same \x) = 
\y) as in Eq. but for 

T=idiag([l,l,0]). (71) 

We find that the CQ/QC discords are Da = Db = 
(3 — -\/3)/64 as in the former example, the CC discord 
is given by Ds = 0.02322 • • • , the nonadaptive upper 
bound is 1)^™''^ = (28 - IIV^) /384 = 0.02330 • • • , and 
the adaptive upper bound is Z?^^"''' = [57 — ll\/3 — 

^^6(62 3v^)] /768 = 0.02326 • • • . The iteration proce- 
dure converges to Ds but not that fast as in the former 
example (see Table I for details). 

Example 3. — Now, let us analyze a state with \x) ^ \y) 
as defined by 

|a:) = i[l,2,3]^, |2;)-^|x), T-idiag([l,2,3]). (72) 
6 7 8 

Analytical formulas for the discords and their upper 
bounds are quite lengthy for this state, so we give only 
their approximate numerical values. The QC discord is 
Db « 0.0259, the CQ discord is Da ~ 0.0262, the CC 
discord is Ds ~ 0.0280, the nonadaptive upper bound 
is d'^s^^^ ~ 0.0284, and the adaptive upper bound is 
d'-s''^'^ « 0.0281. The adaptive upper bounds d'^"^"'^ 
after the nth-iteration are shown in Table I. In conclu- 
sion, we have the following inequalities 

Db<Da<Ds^ i^^""'^^)-O(10-i°) < i?^""'^) < i?^""''^ 

(73) 

All the above examples show how fast the iterations ap- 
proach the correct values of the CC discord. Now, we 
give a counterexample: 

Example 4- — The iteration procedure fails, e.g., for 
the state, given by Eq. (gni for 1/2 < p < 3/5 (see Fig. 
2), as A„ = Aq > for n = 1, 2, .... In general, this can 
be explained as follows: 



Criterion. — If, for a given two-qubit state, the nth iter- 
ation of the adaptive upper bound D^g^^ (in particular 
for n = 0) differs from the CC discord Ds, and |fci-"^) = 
|fj"^) and = 1^1"^) then the iteration procedure 

does not converge to Z^s as = ^ Ds 

for k = 1,2,.... 

Finally, we note that this iteration procedure can be 
improved by replacing D^^"^ by the optimized D^g^^^ 
as described in the preceding subsection. 



VII. CONCLUSIONS 

We have shown that the geometric measures of quan- 
tum correlations, i.e., the CC, CQ and QC discords, are 
equal to the minimal purity deficit under specific von 
Neumann's measurements compatible with the CC, CQ 
and QC classes of states, respectively. This allowed us 
to quickly reproduce known results in the case of qubits 
and also to give some strong arguments that, the CC 
discord may not, in general, be described analytically 
even for a two-qubit state. The best general two-qubit 
formula, given by Eqs. (|Bl[) - (|B6p . still requires minimal- 
ization over two variables. This is in contrast to the 
CQ/QC discords for which analytical two-qubit formulas 
are available. Therefore, we focused on analytical ap- 
proximations of the CC discord. We proposed nonadap- 
tive (i.e., simple product) and adaptive upper bounds for 
the CC discord and applied them for two-qubit states. 
We showed that they are tight and faithful, so they can 
be used as independent tests of nonclassical quantum cor- 
relations. The adaptive upper bound corresponds to an 
optimal measurement on one of the parties conditioned 
an optimal measurement on the other party. We also de- 
scribed a method of improving the adaptive upper bound 
by nonoptimal single-party measurements. This refined 
bound gives exact values of the CC discord for (probably) 
all classes of states for which there are known analytical 
expressions. For randomly generated states, the bound 
usually differs from the CC discord by the order 10"'' 
or 10~^. Moreover, we described an iterative procedure 
for the adaptive upper bound, which usually quickly con- 
verges to the CC discord. We believe that this estimation 
of the symmetric discord will play a role in analyzing the 
cases when all the subsystems of a given quantum system 
interact with the environment on equal footing. For those 
cases it will probably be more adequate than asymmetric 
discord that is based on system-apparatus picture. 

Acknowledgments. This work was supported by 
the Polish National Science Centre under grants 
2011/03/B/ST2/01903 and 2011/02/A/ST2/00305. 



[1] R. Horodecki, P. Horodecki, M. Horodecki, and K. [2] H. Ollivier and W. Zurek, Phys. Rev. Lett. 88, 017901 
Horodecki, Rev. Mod. Phys. 81, 865 (2009). 



10 



(2001). 

[3] L. Henderson and V. Vedral, J. Phys. A 34, 6899 (2001). 

[4] K. Modi, A. Brodutch, H. Cable, T. Paterek, and V. 
Vedral, e-print arXiv: 11 12.6238 (2011). 

[5] J. Oppenheim, M. Horodecki, P. Horodecki, and R. 
Horodecki, Phys. Rev. Lett. 89, 180402 (2002). 

[6] M. Horodecki, P. Horodecki, R. Horodecki, J. Oppen- 
heim, A. Sen De, U. Sen, and B. Synak, Phys. Rev. A 
71, 062307 (2005). 

[7] K. Modi, T. Paterek, W. Son, V. Vedral, and M. 
Williamson, Phys. Rev. Lett. 104, 080501 (2010). 

[8] B. Dakic, V. Vedral, and C. Brukner, Phys. Rev. Lett. 
105, 190502 (2010). 

[9] M. Shi, F. Jiang, and J. Du, e-print llarXiv: 1107.29581 

(2011) . 

[10] S. Luo and S. Fu, Phys. Rev. A 82, 034302 (2010). 
[11] D. Girolami and G. Adesso, Phys. Rev. Lett 108, 150403 

(2012) . 

[12] A. Miranowicz and A. Grudka, J. Opt. B 6, 542 (2004). 



Appendix A: Proofs of observations and lemmas of 
Section II 

Proof of Observation 1. — Since Mx{p) C fix, by def- 
inition, the right-hand-side of Eq. ^ is not greater than 
that of Eq. ([3]). However, one can show that the opposite 
ordering of these expressions can occur hence proving Ob- 
servation 1. Indeed, choose any measurement Aix com- 
muting with (7^. Then the difference of Eqs. ^ and ^ 

is equal to S{p\\a*x) - S{p\\Mx{p)) = S{Mx{p)\W*x) > 
0, where the first equality is due to Lemma 1 below. 
Therefore S{Mx{p)\\(^*x) = = M.x{p)^ where 

the measurement is the optimal one. Lemma 1 further 
furnishes the equivalence between the deficit Ax [p) and 
D^{p) of Eq. dSl). 

Proof of Lemma 1. — Consider any von Neumann's 
measurement M.{-) = 'Yl,iPi{')Pi fo^' orthogonal pro- 
jectors {Pi}, J2iPi — ^- Since any function of 
a Hermitian operator commutes with the operator 
itself, one obtains f{M{G)) = MfiMiG)) and, 
consequently, ti[Ff{M{G))] = ti[FM{fiM{G))] = 

ty[FJ:^P^f{M{G))P,] = tl■[J:,P^FP^f{M{G))] = 

tr[X(F)/(7W(G'))] for any Hermitian F and G. 

Proof of Observation 2. — To prove this Observation 
for the symmetric discord Ds, we consider 

Dsip) = trip^) + min [tT{a^) - 2ti{pa)] (Al) 

optimized over all the CC states a with eigenvectors 
formed by two orthonormal bases Ba ^ Bb ■= {\si)\fj)} 
and some eigenvalues {pij} = ps- We may rewrite 
Eq. (jAll) explicitly as 

Dsip) = tr(p') - ^ min [2tr(pa) - tiia^)]. (A2) 

Ps,Ba<>(jBb 

Let B'^ (E) B*g be an optimal basis defining naturally the 
von Neumann measurement M.\® M.*^ = M.*s- The 
variational state defined in this basis of course satisfies 
M*s{a) = a and so (by Lemma 1) tr(pa) = iT[M*s{p)(T]. 



Denoting by qs the diagonal of the state J^*s{p) (which 
is already of the CC type), we obtain Ds{p) = tr(p^) — 
{2qsPs ~ Ps)' which yields the optimal spectrum 
Ps — Qs- This concludes the proof that the optimal state 
a* in Eq. (jAip satisfies a* — M*g[p). This combined 
with Lemma 1 implies both Eqs. ([S]) and ([5]). 

Consider now the asymmetric discord Da, which is 
given by 

Da{p) = min min [iy{p^) - 2iv{pa) + tr(cr2)],(A3) 

{pi,(Ti} eA = {ei} 

where a = ^iPi\ei){ei\ ® (Ji- Let B'^ = {|e,*)} be an 
optimal basis in Eq. (jA3|) defining now the von Neu- 
mann measurement M*a and the partially optimized class 
of states a' = J2iPi\K){K \ ® which clearly satisfies 
[M*A (E) Ib]{o-') ^ o'. By Lemma 1, 

Da{p) = tr(p2) + min Y}p'M{<f ) - ^p^q^tMp'dl 

(A4) 

where the parameters come from a new state 

p' EE [M*A^lB]ip) = ^/'^ 
i 

For all measurements Mi leaving p'^-s invariant we have 
tr[((T,')2] ^ tr[Mi{a'^f] + Si (with 6^ > 0), since von 
Neumann's measurements do not increase purity. Using 
Lemma 1 again, we therefore have 

min ^[p?tr((a02) - 2p,qM'^'^p'^)] 
- mill Yy>^^^^<f) + ^» - 2Kg,tr(X,(a0pD] 
> min J2[pMM^icT'^f) - 2p,qMM^ia',)p',)] 

= , - .-M, Y.^pM^^'^)^) - 2p^q^Halp'^] 
= min[tr((a')')-2tr(aV')], (A6) 
where 

^'-E^'«l<)(<l®^'(^^)- (A7) 

i 

Since a' and p' commute having product eigenvectors 
(which, however, do not form a product of the two eigen- 
bases in general), we are left only with the final problem 
of finding optimal eigenvalues of a' . In analogy to the 
solution of Eq. (IA2[) . one concludes immediately that 
the optimal spectrum is the same as that of p' . Thus, 
the optimal CQ state must be equal to p', which is just 
the original p subjected to some specific measurement 
M*A ® Ib- This concludes the proof for the CQ-type 
discord and thus, finally, the proof of Observation 2. 

Proof of Lemma 2. — The property (i) is immediate 
since the set Q^s of all the CC correlated states is a sub- 
set of the sets ^a and VIb of the CQ/QC correlated states 
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as shown intuitively in Fig. 1. Property (ii) follows from 
the fact that for Da[p) — the optimal measurement 
M*g providing Db already reduces the state p to a CC 
state. This means that combining it with the measure- 
ment AAa commuting with the left reduction of the state 
yields the same value Db- This value based on product 
measurement, corresponds by definition Eq. ^ to some 
upper bound of the CC discord Ds- On the other hand, 
Db is by (i) a lower bound on the CC discord. Thus 
Ds = Db- Property (iii) is implied by (ii). 



Appendix B: Numerical calculation of the CC 
discord 

For completeness, we give an explicit formula for nu- 
merical calculation of the CC discord for an arbitrary 
two-qubit state p as follows 



Ds{p) = I|p||^-max||cr5(6',( 



^sir, (Bi) 



where as [and thus cr*s — /(ja^s), I2/5), 7^)] can be ex- 
pressed solely in terms of the versor 



\xs) 



= [sin 6 cos (j), sin 9 sin 0, cos 9] ^ 



V{xs\xs) 
as follows 

4\\a*s\\' = l + {x*s\x*s) + {y*s\y*s) + {x*s\n.\y*s} 
= l + {x*s\x) + {y*s\y) + {x*s\T\y*s) 
= 1 + max ((£s|x)2 + {ys\y)^ + {xs\T\ys)^) 



1 + max[Xy{xs) + {xslx^] 
1 + niax[A:r(ys) + {ys\y)^], 



(B2) 



where the quantity Xy(xs) [Xx{ys)] is the maximal 
eigenvalue of the rank- two matrix T\ys){ys\T'^ + 
iT^\xs){xs\T+\y){y\)- So, explicitly. 



and 



where 



\y{xs) - h+ +^{xs\T\y)^ + hl, 
X.iys) - 9+ +^{x\T\ysY+gl, 



h±^ \{{y\y)±{xs\TT'r\xs)), 
g± = \{{x\x) ± {ys\TTT\ys)l 



(B3) 



(B4) 



(B5) 
(B6) 



in agreement with Ref. [9*]. 



Appendix C: Calculation of discords for mixtures of 
|00> and Bell's states 

Here, we give more details of our calculation of the 
CC and CQ/QC discords for the states p{p,(f>) defined 
by Eq. (|46t. 

The correlation matrix T of Bloch's representation for 
the state p{p, (p) reads as 



T 



p cos (p — p sm ( 
p sin (f) p cos (f> 







l-2p 



(CI) 



and the local Bloch's vectors are \x) — \y) — [0, 0, 1— p]'^. 
First, we note that 



K, = Ky = \x){x\+TT^ = 



p2 

p2 

q 



(C2) 



where g = (1 — 2p)^ + (1 — p)^- Since p^ <q\s fulfilled if 
p £ [0, i], so we have to analyze two solutions for p < i 
and p> \ . Thus, the CQ and QC discords are 



Da = Db =- [tr(if:r) - maxeig(i^':r) 



= - [7p2_6p+2-max(g,p2)] (C3) 
= - min(p^, — 3p + 1). 

We can also calculate the CC discord as follows. The 
norm 1 1 175 1 1 ^ is given by 

AWasf - 1 = max {{xs\xY + {ys\yY + {xs\T\ysf) 
xs-ys 



■ max((l|a;) 



(iirii)^ 



{i\xf + {i\yf + {mi?) 

:max [p^2(l -p)2 + (1 - 2pf] . 



(C4) 



2p{p — 1) + 1, since {x\x) = {y\y) — 
2 = (1 - 2pf + 2p2. Thus, finally. 



Moreover, ||p|| = 
(1 - pf and ||T| 
we obtain Ds = \\p\\'^ — Ijfsll^ given by Eq. (ITH) . The 
geometric discords Da ~ Db and Ds for this state are 
plotted in Figs. 2 and 3. 



